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Abstract 

We study the universal properties of the phase diagram of QCD near the critical 
point using the exact renormalization group. For two-flavour QCD and zero quark 
masses we derive the universal equation of state in the vicinity of the tricritical 
point. For non-zero quark masses we explain how the universal equation of state of 
the Ising universality class can be used in order to describe the physical behaviour 
near the line of critical points. The effective exponents that parametrize the growth 
of physical quantities, such as the correlation length, are given by combinations of 
the critical exponents of the Ising class that depend on the path along which the 
critical point is approached. In general the critical region, in which such quantities 
become large, is smaller than naively expected. 



1 Introduction 



The most prominent feature of the phase diagram of QCD at non-zero temperature and 
baryonic density is the critical point that marks the end of the line of first-order phase 
transitions. (For a review see ref. [1].) Its exact location and the size of the critical region 
around it determine its relevance for the heavy-ion collision experiments at RHIC and 
LHC. 

There exist several analytical studies of the phase diagram near the critical point 
[2]-[12]. Lattice simulations are a source of precise information for the details of the 
phase diagram. Recent studies [13, 14, 15] have overcome the difficulties associated with 
simulating systems with non-zero chemical potential, even though the continuum limit or 
realistic quark masses have not been reached yet. 

The critical equation of state encodes all the physical information for a system near a 
critical point in the presence of an external source. In the case of QCD, the parameters 
that can be adjusted in order to approach the critical point are the temperature T and 
baryonic chemical potential [i. For two-flavour QCD, the role of the external source j is 
played by the current quark mass that breaks explicitly the chiral symmetry of the QCD 
Lagrangian. If isospin-breaking effects are neglected, the source is a (linear) function of 
the common mass of the light quarks. 

An analytical description of the critical region requires the use of the renormalization 
group. The study of first-order phase transitions is easier within the Wilsonian (or exact) 
formulation [16]. We shall use the formalism of the effective average action T k [17, 18], 
that results by integrating out (quantum or thermal) fluctuations above a certain cutoff 
k. The dependence of T k on k is given by an exact flow equation [19, 20]. When k is 
lowered to zero becomes the effective action. The effective average action can by used 
in order to desribe very efficiently the universal and non-universal aspects of second-order 
phase transitions, as well as strong or weak first-order phase transitions [21]. 

In order to discuss the phase diagram we have to employ an appropriate order pa- 
rameter. The usual choice is the quark-antiquark condensate that is associated with the 
spontaneous breaking of the chiral symmetry of the QCD Lagrangian. For two flavours 
an equivalent description uses as effective degress of freedom four mesonic scalar fields 
(the a-field and the three pions 7Tj), arranged in a 2 x 2 matrix 

$ = I(a + m-f), (1) 

where Tj (i — 1, 2, 3) denote the Pauli matrices. The various interactions must be invariant 
under a global SU(2)l x SU(2)r symmetry acting on $. If the a-field (that corresponds 
to a condensate uu + dd) develops an expectation value the symmetry is broken down 
to SU(2)l+r- The explicit breaking of the chiral symmetry through the current quark 
masses can be incorporated as well. It corresponds to the interaction with an external 
source through a term — ja in the Lagrangian. 

The full QCD dynamics cannot be described by an effective Lagrangian involving only 
mesonic fields. More complicated effective descriptions, such as the linear quark-meson 
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model, have been used for the discussion of the phase diagram for low chemical potential 
[22, 23, 24]. In ref. [25] it was shown that the expected structure of the phase diagram for 
two-flavour QCD can be reproduced within this model. In particular, in the absence of an 
external source (zero current quark mass), the phase diagram contains a line of first-order 
phase transitions and a line of second-order ones, that meet at a tricritical point. The 
second-order phase transitions belong to the 0(4) universality class, while the tricritical 
point is associated with the Gaussian fixed point and is described by mean-field theory. 
In the presence of a source the first-order line ends at a critical point. The second order 
phase-transition at this point belongs to the Ising universality class. The reason is that 
all degrees of freedom other than the a-field remain massive and decouple at low energy 
scales. For varying external source we obtain a line of second-order phase transitions, that 
starts at the tricritical point and marks the boundary of a surface of first-order transitions. 

The universal properties associated with the lines of second-order phase transitions 
in the phase diagram do not depend on the effective model employed for the description 
of the physical system (the mesonic model described above, a discretized version of the 
QCD Langrangian employed in lattice simulations, etc). The crucial ingredient is the 
inclusion of the degrees of freedom that become critical (massless) at a transition, and 
the symmetries that characterize the effective interactions [26]. For the same reason the 
inclusion of the more massive quarks is not expected to modify the universal behaviour. 
It may only shift the location of the various lines on the phase diagram. 

The purpose of this paper is to derive the universal equation of state of QCD in 
the region of the phase diagram near the line of critical points. We concentrate on the 
dependence on the quantities that can be varied In the following section we desribe a 
simplified model that has a phase diagram similar to that of QCD at relatively small 
values of the chemical potential. In section 3 we discuss the phase diagram in the context 
of mean-field theory. In section 4 we describe the renormalization-group approach. We 
also determine the initial conditions for the renormalization-group flow in the context of 
the quark-meson model. In section 5 we derive the equation of state in the absence of 
an external source. In section 6 we derive the equation of state for a non-zero external 
source in the vicinity of the line of critical points. In the final section 7 we discuss the 
implications of our results for the possibility of observing the universal critical behaviour 
of QCD in heavy-ion experiments. 

2 The model 

As we are interested only in the universal characteristics near the second-order phase 
transitions, we need only discuss a simplified model that preserves the basic structure of 
the phase diagram of two-flavour QCD. We consider a Lagrangian in three dimensions 
that contains standard kinetic terms for the fields a, 7Tj (i = 1, 2, 3), and the potential 

V{p)=m 2 p+ l -\p 2 + l -gp\ (2) 
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with p = Tr = (a 2 + 7T i 7T t )/2. The couplings m 2 , A, g are considered functions of 

the temperature T and the chemical potential /i. We assume that g is always positive, so 
that the potential is bounded from below. 

In ref. [25] it was shown that such a potential emerges in the quark-meson model after 
the integration of the temperature fluctuations with non-zero Matsubara frequences. The 
fermions are completely integrated out as their lowest Matsubara frequency is non-zero. 
The only mode that survives at energy scales below the temperature is the zero-mode of 
the mesonic field. Thus the theory becomes effectively three-dimensional. The effective 
potential around the origin has a structure similar to that of eq. (2). In section 4 we 
calculate the effective potential in the context of an effective model of QCD, the linear 
quark-meson model. 

The various fixed points that determine the structure of the phase diagram appear 
during the integration of the low-energy fluctuations of the zero mode. Theories for which 
the renormalization-group flow is dominated by an infrared fixed point become insensitive 
to the details of the ultraviolet formulation. For this reason the simplified model we are 
considering is sufficient for the discussion of the universal properties of the phase diagram 
of QCD near the lines of second-order phase transitions. 

The effect of a current quark mass can be taken into account by adding a term —jcr, 
with j a (linear) function of the mass. 

3 Mean-field theory 

It is instructive to discuss the phase diagram of our simplified model neglecting the field 
fluctuations. We allow for variations of m 2 and A, while we assume that the coupling g 
remains constant. This is a good approximation of the behaviour of the potential studied 
in ref. [25]. 

We begin by studying the phase diagram in the absence of an external source (j = 0). 
Let us consider first the case A > 0. For m 2 > the minimum of the potential is located 
at a = iii = and the system is invariant under an 0(4) symmetry. For m 2 < the 
minimum moves away from the origin. Without loss of generality we take it along the 
a axis. The symmetry is broken down to 0(3). The 7r's, which play the role of the 
Goldstone fields, become massless at the minimum. If m 2 (T,fi) has a zero at a certain 
value T = T cr , while X(T cr ,fi) > 0, the system undergoes a second-order phase transition 
at this point. The minimum of the potential behaves as a ~ |T — T cr \^l 2 slightly below 
the critical temperature. The critical exponent (3 takes the mean-field value (3 = 1/2. 

Let us consider now the case A < 0. It is easy to check that, as m 2 increases from 
negative to positive values (through an increase of T for example), the system undergoes 
a first-order phase transition. For |A| approaching zero the phase transition becomes 
progressively weaker: The discontinuity in the order parameter (the value of a at the 
minimum) approaches zero. For A = the phase transition becomes second order. The 
minimum of the potential behaves as a ~ |T — T^l 1 / 4 . The critical exponent (3 takes the 
value (3 = 1/4 for this particular point. 
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Figure 1: The phase diagram. 



Let us assume now for simplicity that A is a decreasing function of fi only, and has a 
zero at fi = /i*. If we consider the phase transitions for increasing T and fixed \x we find 
a line of second-order phase transitions for //<//*, and a line of first-order transitions for 
ji > /i*. The two lines meet at the special point (T*, /i*), where T* = T cr for /i = //*. This 
point is characterized as a tri critical point. (For a review of the theory of tricritical points 
see ref. [31].) In the general case A will be a decreasing function of a linear combination 
of fj, and T. The structure of the phase diagram remains the same. Simply there is a 
linear combination of T and /i that generates displacements along the lines of first- and 
second-order phase transitions near the tricritical point, and a different one that moves 
the system through the phase transition. 

If a source term — ja is added to the potential of eq. (2) the 0(4) symmetry is explicitly 
broken. The phase diagram is modified significantly even for small j. The second-order 
phase transitions, observed for A > 0, disappear. The reason is that the minimum of the 
potential is always at a value a ^ that moves close to zero for increasing T. For small 
j the mass term at the minimum approaches zero at a value of T near what we defined 
as T cr for j = 0. However, no genuine phase transition appears. Instead we observe an 
analytical crossover. 

The line of first-order transitions persists for j ^ 0. However, at the critical temper- 
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ature the minimum jumps discontinuously between two non-zero values of a and never 
becomes zero. The line ends at a new special point, whose nature can be examined by 
considering the a-derivative of the potential V(p) of eq. (2). For a first-order phase 
transition to occur, dV/ da must become equal to j for three non-zero values of a. This 
requires d 2 V/da 2 = at two values of a. The critical point corresponds to the situation 
that all these values merge to one point. At the minimum a* of the potential V(a;j) = 
V(a,7Ti = 0) — ja at the critical point, we have: dV /da = d 2 V /da 2 = d 3 V/da 3 = 0. It 
can be checked that this requires A < and can be achieved by fine-tuning m 2 and A 
for given j. The minimum cr* is then completely determined. At the critical point we 
expect a second-order phase transition for the deviation of a from a*. As d A V/da 4 ^ at 
the critical point, the critical exponent (3 takes the value (3 = 1/2. The 7r's are massive, 
because m 2 = d 2 V(p)/dn 2 = dV(p)/dp = j /a* at the critical point. 

The phase diagram we discussed is presented in fig. 1. If we take into account the 
effect of fluctuations of the fields, the potential will have a more complicated form than 
the simple assumption of eq. (2). As we shall see the qualitative structure of the phase 
diagram is not modified. However, the nature of the fixed points may differ from the 
predictions of mean-field theory. 



4 The renormalization-group flow 

The most efficient way to study a physical system near a phase transition is through 
the effective potential. In the formulation of the renormalization group which we are 
employing, the dependence of the potential of a three-dimensional theory on a "coarse- 
graining" scale k is described by the equation [19, 21] 

^u k {a) = -3u k + 1(1 + V )au' k + ^ [Z 3 «) + 3/ 3 (u'Ja)} , (3) 
where t = ln(/c/A) and we have defined the dimensionless quantities 

u k = k~ 3 U k , a = k~^z\a. (4) 

Primes denote derivatives with respect to a. The scale-dependent potential U k results 
from the integration of the field fluctuations with characteristic momenta larger than k. 
In the limit k — > it becomes equal to the effective potential. It is possible that the 
three-dimensional theory described by the potential U k has its origin in the dimensional 
reduction of a more complicated, higher-dimensional theory. This situation arises during 
the study of the phase diagram of QCD at non-zero temperature T and baryonic chemical 
potential p (see below). 

The above equation can be derived from an exact flow equation if only two terms 
are retained in the effective action: the potential and a kinetic term that includes a 
$- independent wavefunction renormalization Z k [21]. The anomalous dimension rj = 
—d(lnZ k )/dt describes the scale dependence of Z k , and can be determined starting from 
the exact flow equation. It becomes constant when the flow approaches a fixed point [21]. 
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The "threshold" function Zjj is a particular case of 

l» = l^-^-iT J d d q ^(P k + wkr (n+1) . (5) 

The inverse propagator for massless fields has been modified so at to forbid the propa- 
gation of modes with momenta below the "coarse-graining" scale k. Several forms of the 
modified propagator are possible, which correspond to sharper or smoother cutoffs of the 
low-momentum modes [21]. 

The continuous reduction of k, from the ultraviolet scale A to zero, permits the gradual 
integration of fluctuations with different characteristic momenta. The function Iq(w) falls 
off for large values of w following a power law. As a result it introduces a threshold 
behaviour for the contributions of massive modes to the evolution equation. When the 
running squared mass of a massive mode becomes much larger than the scale k 2 , the mode 
decouples and the respective contribution vanishes. In eq. (3) we distinguish contributions 
from two different types of fields: the radial mode (the a-field) and the three Goldstone 
modes (the pions). Their masses are expressed through the derivatives of the effective 
potential. 

The initial condition that is needed for the solution of the evolution equation (3) is 
given by eq. (2). The potential at the ultraviolet scale A is identified with the classical 
(or bare) potential: U\ = V. The integration of field fluctuation with momenta between 
A and zero leads the determination of the effective potential: U e ff = Uo = U. 

It is possible that the theory at k — A results from the integration of some complicated 
ultraviolet theory. This situation arises during the study of the phase diagram of QCD 
at non-zero temperature T and baryonic chemical potential p [1, 23, 24, 25]. An effective 
theory that captures many of the features of QCD is the quark-meson model [22]. It 
employs mesons and quarks as physical degrees of freedom. The QCD phase transition is 
associated with the breaking of the chiral symmetry through a non-zero expectation value 
for the mesonic a-field. The connection with the model we are studying is made through 
the phenomenon of dimensional reduction at non-zero temperature. At energy scales 
below the temperature a four-dimensional theory can be described in terms of effective 
three-dimensional degrees of freedom. In a Fourier decomposition of a scalar field only 
the term with zero Matsubara frequency (the zero mode) survives. The other modes 
decouple, as they get a mass proportional to the temperature. The fermions do not have 
a zero mode and they decouple completely. 

In the context of the quark-meson model with two flavours an analytical expression 
has been derived for the form of the potential at a scale k = A = T [25]. It is 

U' A (p;T,p)~\ 0R 

where 

h 4 I & \ 

I F (p; T, p) = ~N c — 2 p In + V F (p; T, p) + V F (p; T, -p) , (7) 

71 V + 2A 2 / 



P ~ POR ~ 



An 2 ' 



e~ 2 



7T 

~3~ 



+ I F (p;T, p), 



(6) 
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V F (p;T,p) = -4N C T J -0-\n(l+exp 

— oo 



~^^q 2 + h 2 p/2-p 



(8) 



Here pqr is a renormalized expectation value for the mesonic field at zero temperature 
and chemical potential, while Ao_r is the renormalized mesonic quartic coupling. The term 
~ T 2 in eq. (6) results from the integration of the scalar thermal fluctuations with non- 
zero Matsubara frequences. The effect of the fermionic quantum and thermal fluctuations 
is incorporated in the quantity Ip(p;T,p). Even at zero temperature, the fermionic 
quantum fluctuations give a contribution (the first term in the r.h.s. of eq. (7)) that has 
to be taken into account when defining the vacuum of the zero-temperature theory. The 
quantity Vf(p;T, p) is the free energy of a fermionic gas at non-zero temperature and 
chemical potential. We have taken into account only two quark flavours (corresponding 
to the light u- and o?-quarks), as well as their antiparticles through V F (p;T, —p). The 
parameter N c = 3 stands for the number of colours. The scale A corresponds to the 
energy scale at which the quark- meson model becomes an effective description of QCD. 
It is close to 650 GeV [22], while T ~ 200 MeV. The constituent quark mass is given by 
h 2 p/2. 

We emphasize that the above expressions have been derived for small A _r and h. If we 
would like to make contact with realistic mesonic physics we must reproduce the correct 
pion decay constant /„., the a-field mass and the constituent quark mass at zero temper- 
ature and chemical potential. If we ignore the logarithmic contribution in Ip(p;T, p), we 
have to set \ / 2p 0R ~ f n 
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87 MeV, m a ~ V^orPor - 600 MeV, M q = ^h 2 p m /2 ~ 300 
MeV. This leads to A, h — O(10). For this reason the essentially perturbative expressions 
(6)-(8) must be used only for qualitative discussions of the behaviour of realistic QCD. 
Whenever we make use of these expressions in the following we assume that \or and h 
are small, consistently with their derivation. 

The quantity I F (p] T, p) can be expanded in powers of p as 



I F (p; T,p) 



h 2 N c 
12 



T 



7T 



-p 



h 4 N c 
16tt 2 



-1 + 2 7e + In 



4T 2 



+ 2Li( 1 >°) (0,-exp ^)+2Li( 1 '°) (o, — exp — 



where 7 e ~ 0.5772 is the Euler-Mascheroni constant and Li^ ll ' l2 \n,x) denotes the l±-th 
and /2-th partial derivative of the polylogarithmic function Li(n, x) with respect to n and 
x respectively. For p = we have Li^'°^(0, — 1) = — ln(7r/2)/2 and we recover the known 
expressions of ref. [27]. 

Making use of eqs. (6), (9) we can write the potential Ua in the form of eq. (2) with 1 



m 



— AoR P0R + 



A, 



OR 



1 - 



7T3/2 



+ 



h 2 N r 



12 



T 2 + 



h 2 N, 



4tt 2 



c 2 
p 



(10) 



1 The quartic coupling of the effective three-dimensional theory is given by AT. We do not refer to 
this coupling, as its additional temperature dependence is eliminated trivially when we discuss the phase 
diagram in terms of the observables of the four-dimensional theory 
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Figure 2: The evolution of the rescaled potential for m 2 = 0.0213, A = —0.335, g — 1. 



A 



Aqr + 



h 4 N c 
16tt 2 



-1 + 2 7e + In 



4T2 



+ 2 Li™ (0,-exp(|)) 

+ 2 Li™ (0,-exp (-£)) 



.(11) 



The function Li*- 1 ' - 1 (0, — exp(x)) + Li*- 1,0 - 1 (0, — exp(— x)) is monotonically decreasing, and 
it takes the value — ln(7r/2) for x — 0. These expressions demonstrate how the phase 
diagram that we discussed in the previous section can emerge in the context of QCD. 
Even if the mass term is negative at T — fj, — 0, it can become positive for sufficiently 
large T or fi. For constant T, the quartic coupling is a decreasing function of fi. For small 
h the tri critical point is expected to appear for values of T and \i such that A ~ 0. 



5 The equation of state near the tricritical point 

5.1 The crossover behaviour 

The phase structure in the absence of an external source (j = 0) can be determined in 
complete analogy to ref. [25]. There is a line of second-order phase transitions which can 
be approached for fixed positive values of A, g in eq. (2) by fine-tuning the mass term. 
During the renormalization-group flow the potential approaches a form characteristic of 
the fixed point of the 0(4) universality class. Subsequently the system deviates towards 
the broken or symmetric phase. The influence of the fixed point results in a behaviour that 
can be characterized by universal quantities, such as critical exponents. For sufficiently 
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Figure 3: The effective critical exponents v and (3. 

negative values of A the phase transitions become first-order. Again, an appropriate fine- 
tuning of the mass term leads to an effective potential with two coexisting phases. No 
fixed point is approached during the flow. 

The lines of first- and second-order phase transitions meet at a point characterized by 
specific values of m 2 and A. The evolution of the potential in the proximity of this point 
is depicted in fig. 2. The initial condition is given by eq. (2). For the plot we integrated 
eq. (3) with rj = (and, therefore, Z k = 1). We have used a smooth infrared cutoff in the 
propagator P k of eq. (5), similarly to ref. [25]. The anomalous dimension is very small 
[r] ~ 0.03) in this model and setting it equal to zero gives a good approximation to the 
exact solution. In fig. 2 we plot the quantity u' k /a as a function of a 2 /2 for Hi = and 
decreasing values of the scale k. In an alternative notation, that makes the 0(4) symmetry 
apparent, the figure depicts [dU k j 'dp)/k 2 as a function of p — p/k, with p defined below 
eq. (2). 

Two quantities must be fine-tuned for the type of evolution shown in fig. 2. We fixed 
g — 1 and varied A. (All quantities are given in units of the ultraviolet scale A.) For every 
value of A we fine-tuned m 2 so as to be very close to the phase transition. For A = —0.335, 
m 2 cr ~ 0.0213 (the latter quantity is determined with a precision of 15 significant figures) 
we are in the immediate vicinity of the tricritical point and very close to a second-order 
phase transition. 

Initially the potential is very flat in the range of a we are considering, so that its 
derivative is almost zero. The system is very close to the Gaussian fixed point. As k 
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Figure 4: The evolution of the rescaled potential for m 2 = 0.0217, A = —0.338, g — 1. 

becomes smaller another fixed point is approached, characteristic of the 0(4) universality 
class. Eventually the system moves towards the symmetric phase (u' k /a becomes positive 
at the origin). In the limit k — > we obtain the effective potential, from which we can 
derive renormalized quantities such as the location of the vacuum and the field mass. 

For systems close to a second-order phase transition, we can parametrize physical 
quantities in terms of critical exponents. In our example, we can keep A and g fixed with 
the values A = —0.335, g — 1 and vary the mass term according to m\ = m 2 cr + 5m 2 . For 
5m 2 > the renormalized mass term at the origin can be parametrized as m\ ~ (5m 2 ) 2u . 
For 5m 2 < the minimum of the potential is located at a non-zero field value that can 
be parametrized as cr 0R ~ |5m 2 |^. (The quantity 5m 2 replaces the quantity T — T cr 
that appears when the phase transition is discussed in the context of the fundamental 
four-dimensional theory at non-zero temperature.) 

The effective critical exponents is, (3 are depicted in fig. 3 as a function of 5m 2 . We 
observe that for relatively large 5m 2 they take the values (3 = 0.25, v = 0.5. These are 
the mean-field predictions that we derived in section 2. They are the expected values of 
these exponents, as predicted by the general theory of tricritical points [31]. For large 
5m 2 the potential never approaches the 0(4) fixed point. It stays close to the Gaussian 
one, before moving directly towards the broken or the symmetric phase. For smaller 5m 2 
system feels the attraction of the 0(4) fixed point. For 5m 2 — > the potential spends a 
large part of its evolution near this fixed point. As a result it loses memory of its initial 
form near the Gaussian fixed point. The critical exponents take the values characteristic 
of the 0(4) universality class (3 = 0.41, v = 0.82. 2 The curves of fig. 3 are typical 

2 The most accurate known values for these exponents are close to (3 = 0.39, v = 0.75 [37]. The values 
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examples of crossover curves. They describe the variation of universal quantities, such as 
the critical exponents, as the relative influence of two fixed points changes. 

In fig. 4 we display the evolution of u' k /a for a theory with m 2 ~ 0.0217, A = —0.338, 
v — 1. The slight decrease of the value of A relative to fig. 3 modifies drastically 
the evolution. Initially the potential stays very close to the Gaussian fixed point, but 
eventually moves away from it. Its final form is characteristic of a first-order phase 
transition as was discussed in detail in ref. [25]. This transition is very weak, as the 
discontinuity in the order parameter is several orders of magnitude smaller than A. 

The tricritical point is located on the plane (m 2 , A) in the region between the points 
(0.0213, -0.335) and (0.0217, -0.338) for g = 1. For the critical values m 2 and A t the 
solution of the evolution equation stays close to the Gaussian fixed point for an infinitely 
long "time" t. As a result the critical behaviour is determined completely by mean-field 
theory. For example, the critical exponent (3 stays equal to 0.25 arbitrarily close to the 
phase transition. The relative influence of the two fixed points depends on the magnitude 
of two parameters that can be taken as 5m 2 and 5\ = A — X t . For fixed 5X and Sm 2 — > 
the 0(4) fixed point dominates, while for fixed Sm 2 and — > only the Gaussian fixed 
point is approached. All this is consistent with the general theory of tricritical points [31]. 

5.2 The universal equation of state 

All the physical information near a second-order phase transition can be encoded in the 
equation of state. This describes the relation between the order parameter (in our case 
the field expectation value), the parameters that control the distance from the phase 
transition (such as 5m 2 ) and the external source. In the case of one relevant parameter 
the equation of state can be cast in the form [28, 3] 

with U(a) the effective potential. For sufficiently small 5m 2 , so that the 0(4) fixed point 
is approached, the function f(x) takes a universal form fo(4)(x) that has been computed 
in ref. [24]. The critical exponents 5, (3, as well as the other critical exponents and 
amplitudes encoded in f(x), take values characteristic of this universality class: 5 ~ 4.8, 
(3 ~ 0.39. (Our approximate calculation above gave (3 ~ 0.41). 

For a theory that only approaches the Gaussian fixed point the equation of state can be 
inferred more easily. As we saw earlier, because of the smallness of the various couplings 
near the Gaussian fixed point, the mean- field description is sufficient. For a more detailed 
understanding of the renormalization-group flow, we assume that the potential Uk can be 
approximated by a polynomial, similarly to eq. (2) 

U k (p) = m 2 (k)p + iA(A;)p 2 + ±g(k)p 3 + lg(A;)p 4 + ... (13) 

we are quoting agree with the results of an analysis based on the exact renormalization group within the 
limited truncation we are employing [21]. 
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Starting from the flow equation (3) we can derive evolution equations for the various 
parameters [20] 



dm 2 (k) 
dk 

d\{k) 
dk 




^¥ 2 \¥)~n 
15 A 3 3 (m 2 \ 




(15) 



(14) 



dg{k) 
dk 



(16) 



(The anomalous dimension is zero near the Gaussian fixed point.) The most important 
terms in the above equations are the term in the r.h.s. of the first equation and the 
last term in the r.h.s. of the second equation. These generate contributions to the 
renormalization of m 2 and A that are proportional to the ultraviolet scale A. The last 
term in the r.h.s. of the third equation does not generate significant contributions because 
h = in the bare potential of eq. (2). 

For a potential that is given by eq. (2) for k = A, and a flow that stays in the 
vicinity of the Gaussian fixed point, the integration of the evolution equations is expected 
to renormalize significantly only m 2 and A. The coupling g is expected to receive only 
logarithmic corrections (it corresponds to a marginal operator) cut off by the mass term. 3 
Moreover, these corrections are expected to reduce g(k) and shift the potential closer to 
the Gaussian fixed point [29] . These conclusions are verified by the numerical integration 
of the full evolution equation (3). 

If we neglect the logarithmic corrections a good approximation of the effective potential 

is 



The tricritical point corresponds to the theory with \r — 0, for which the equation of 
state can be written as 



Comparison with eq. (12) indicates that for the tricritical point 5m 2 = 4m 2 R /g, Sj = 4j/g, 
5 — 5, (5 — 1/4, ft r (x) = 1 + x. This is in agreement with the values of the exponents in 
fig. 3 for relatively large 8m 2 . 

For Ar ^ the equation of state can be written in terms of the potential of eq. 
(17) with the identification 5m 2 = m 2 R) 8\ = \r. We expect a crossover away from the 
tricritical point for 5 A > 5m 2 . This is apparent in fig. (3) in the transition of the effective 
exponent f3 from the tri critical-point value f3 = 0.25 to a value close to 0.5. Of course, the 

3 It must be pointed out that no significant infrared effects are expected near the Gaussian fixed point 
because of the decoupling of the various modes at low enough scales. This happens in the symmetric 
phase, around which we expanded the potential in eq. (13). However, if the vacuum is at a non-zero 
value of the field the contributions of the three massless pions (the Goldstone modes) never decouple. 
This leads to interesting effects that have been discussed in ref. [30], for example. 




(17) 
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asymptotic value is determined by the 0(4) fixed point that attracts the flows after they 
leave the Gaussian fixed point. In analogy with the crossover observed in fig. 3 for the 
critical exponents, we expect the critical equation of state to change continuously between 
the form characteristic of the tricritical point and the one of the 0(4) universality class, 
as A is varied. 



6 The equation of state near the critical point 
6.1 The vicinity of the tricritical point 

As we have seen in section 3, for non-zero external source j and negative A it is possible to 
have a second-order phase transition at the end of a line of first-order phase transitions. 
The order parameter is the deviation of the field from a certain non-zero expectation 
value. This type of transition is not associated with the breaking of a symmetry, as the 
original symmetry is explicitly broken by the source. 

For the part of the line of critical points in the vicinity of the tricritical point the 
potential is approximately given by eq. (17) (neglecting logarithmic corrections) and the 
analysis of section 3 based on mean-field theory is sufficient. As we saw there, we must 
define a new effective potential U(a;j) = U(a,7ii = 0) — j*cr whose minimum is located 
away from the origin. The value of the source j* must be defined appropriately. 

In order to derive the equation of state we shift the field: a — > a* + a. The value of 
a* can be determined by requiring that at the critical point the minimum of the potential 
is located at a — 0. As we saw in section 3, at the critical point we have dll /do = 
d 2 U/da 2 = d 3 U/da 3 = at the minimum. The last equality leads to a 2 = — 3A_r/(5#), 
indicating that the existence of a critical point requires \r < 0. 

On a surface of constant \r and cr* (we assume a fixed renormalized g throughout the 
paper) the effective potential for the shifted field can be written as 

^ = ~d~o = l ga * + a * 6rh2 + 6rh2a + t ga * a3 + l 9a * aA + \ 9<T *> ( 19 ) 

where we have redefined the mass term according to m 2 = ml + 5m 2 , with ml = 5gaf/4. 
If we repeat the calculation for a different value of \r, the form of the above equation 
remains the same. The only parameter that changes is a* through its dependence on \ R 
given above. The field a corresponds now to deviations from the new value of cr*. 

It is convenient to define cr* = o"* — e, with being kept fixed and e accounding for 
variations of \r around a constant value. Moreover, the field a can be redefined as the 
deviation from c^o- Then, eq. (19) retains its form with the replacements a* — > — e, 
a — > a + e, and can be written as 

5j = F x (e, 5m 2 ) + (a + e) 3 [j^y 2 + l) + O ((a + e) 4 ) , (20) 

where j = + 5j, = 2gal /3, and 5j = 35j/(5ga 2 ), 5m 2 = 3 5m 2 /(5ga 2 ). We 
assume that |er| takes values in a range of a few |e| <^ a* . The function iq(e, 5m 2 ) 
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is zero for e = 5m 2 = 0. Near the tricritical point it takes the form Fi(e,Sm 2 ) ~ 
-2a 2 e + 4(j, e 2 - 4e 3 + a^5m 2 . 

It must be emphasized that the above expression is a rather complicated way to rewrite 
the potential of eq. (17). As a result, the phase diagram as a function of 5m 2 , 5j and 
e is expected to have the same structure as the one discussed in section 3. This can be 
verified easily. There is a line of critical points parametrized by e for 5m 2 — 8j — 0. A 
particular point can be approached on a surface of constant e. For example, on the surface 
e = the critical point can be approached along the line Sj = i<i(0, 8m 2 ) = a^ 5m 2 . The 
resulting second-order phase transition is parametrized by the critical exponents of mean- 
field theory: j3 = 1/2, 7=1. For e = 0, 5m 2 < the potential has two minima and 
there are two possible vacuum states. A first-order transition from the vicinity of one 
vacuum to the other takes place when 5j is varied through a* 5m 2 . For 5m 2 > the 
transition from positive to negative field expectation values is continuous, but the mass 
term never becomes zero. This situation corresponds to an analytical crossover. The 
qualitative picture remains the same for other values of e. 

The term Fi(e,5m 2 ) in eq. (20) can be absorbed in the source term 5j through a 
redefinition of the potential. However, our choice of 5m 2 , 5j gives the most transparent 
connection with the parameters used for the phase diagram of QCD. The term 5j corre- 
sponds to deviations of the explicit symmetry-breaking term in the bare Lagrangian from 
some constant value. In this sense it corresponds to variations of the current quark mass. 
The terms e, 5m 2 are related to parameters of the theory that do not break explicitly 
the original symmetry. In the QCD case, they are functions of the temperature and the 
baryonic chemical potential. 

Within the quark-meson model, that constitutes an effective description of low-energy 
QCD, the parameters e, 5m 2 can be expressed as functions of 5T and 5fi, where 5T = 
T — T*, 5fi = fi — fi* are the deviations from the values that correspond to the second-order 
phase transition for a given current quark mass. We can integrate eqs. (14)-(16) from 
k = A = T down to k — using (10), (11) as initial conditions for m 2 and A, and taking 
g = 0(h 6 ) ~ 0. 4 The only significant contribution resulting from this integration is a 
term 3A rT 2 /(47t 3 / 2 ) that must be added to the mass term m 2 [25]. We find 



In the following subsection we shall see that it is possible to cast the equation of state in 
a form similar to eq. (20) even away from the tricritical point. An important observation 
is that there are many ways to approach a critical point. The universal properties of 

4 Here h is the Yukawa coupling that determines the constituent quark mass. We assume that it is 
small, consistently with the derivation of eqs. (10), (11). 
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the physical system (such as the values of the critical exponents) are not the same along 
every path. Contrary to the case discussed above, the critical exponents can take effective 
values different from the expected ones. For example, if the critical point is approached 
on the surface Sj = 0, along a line e = e(5m 2 ) such that Fi(e(5m 2 ), 5m 2 ) = 0, we obtain 
a + e(5m 2 ) ~ l^m 2 ) 1 / 2 for the location of the minimum of the potential when 5m 2 < 0. 
The mass term at the vacuum scales ~ |5m 2 |, and the exponent 7 takes the value 7 = 1. 
If the critical point is approached along a line e = e(5m 2 ) such that Fi(e(Sm 2 ), 8m 2 ) 7^ 0, 
we obtain a + e{5m 2 ) ~ l^m 2 ) 1 / 3 for the location of the minimum of the potential when 
5m 2 < 0. The mass term at the vacuum scales ~ \5m 2 \ 2 ^ 3 , and the exponent 7 takes an 
effective value 7 e // = 2/3. In the case of the quark-meson model, it is obvious from the 
form of eqs. (21), (22) that the generic trajectory that passes through the critical point 
has Fi(e(Sm 2 ), 5m 2 ) 7^ with high probability. 

These findings are consistent with the general analysis of critical behaviour in the 
vicinity of tricritical points [31]. According to this analysis one must introduce more than 
one scaling fields [32, 33, 34], and the critical behaviour depends on their relative magni- 
tude. We emphasize, however, that the values that we derived above do not correspond to 
a tricritical theory. The reason is that we are studying the theory near the line of critical 
points. The part of this line that lies in the vicinity of the tricritical point is described 
by mean field theory. As a result we find the mean-field exponents of the critical and not 
the tricritical theory. 

6.2 Away from the tricritical point 

Far from the Gaussian fixed point the theory is renormalized significantly and the potential 
cannot be approximated by a simple polynomial any more. The flow equation (3) describes 
the evolution of the potential as the scale k is lowered from A to zero. The last term in 
the r.h.s. of this equation includes contributions from the radial mode (the a field) 
and the three Goldstone modes (the pions). In the presence of an external source the 
ground state of the system is located at a point where the potential has a non-zero 
derivative and satisfies dll/da = j*. For zero anomalous dimension, the argument of the 
second "threshold" function is expected to become asymptotically u' k /a — > (j*/o"*)//c 2 for 
k — > 0. As a result the pions are expected to decouple and their contribution to the 
evolution to switch off. (A small anomalous dimension does not affect this conclusion.) 
This decoupling behaviour has been demonstrated explicitly in ref . [25] , and we shall not 
repeat the analysis here. 

Deep in the infrared, after the pion decoupling, the flow equation takes the form 

^u k (a) = -3u k + ^(l+v) °< + ^ l 3 o K) ■ (23) 

Apart from the Gaussian fixed point, the only known t-independent solution of this equa- 
tion is the Ising fixed point, for which dui/dt = [36]. This solution has a Z2 symmetry 
a <-> —a. Whether the fixed point will be approached depends on the initial condition for 
the potential. If the potential at some scale k = A after pion decoupling is Z 2 -symmetric, 
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the subsequent flow preserves the symmetry. One needs to fine-tune only one parame- 
ter, the mass term, in order to approach the Ising fixed point. The theory has only one 
relevant parameter. 

However, in the presence of a source the potential U^(a;j) = U^(a, 7ii = 0) — j*cr does 
not possess a Z 2 symmetry and the flow becomes more complicated. If we expand Uk 
after shifting the field by a /c-independent value a — > a* + a, we obtain a polynomial with 
several even and odd powers of a. The crucial question is how many of them need to be 
fine-tuned in order to approach the Ising fixed point. 

In order to answer this question we must study the small perturbations around the 
Ising fixed-point solution ui and examine how many have negative eigenvalues so that they 
grow in the infrared. In linear perturbation theory around the fixed point (uk = uj + Suk) 
the evolution equation becomes 

pu k (a) = -3 5u k + 1(1 + V )a5u' k - JL 6 u' k ' If «') , (24) 

where we have used the property lf(w) = — d[lQ(w)]/dw of the "threshold" functions 
defined in eq. (5). 

There is one even solution 5u±k of the above equation, whose eigenvalue Ai is related 
to the exponent v in the Ising universality class through Ai = \jv. The simplest odd 
perturbation is 5u 2 k = c 2 5", with an eigenvalue A 2 = —(5 — rj)/2. There is one more odd 
perturbation with negative eigenvalue. It is given by Su^k = c^dui/da and its eigenvalue 
is A3 = —(1 + i])/2 [37, 38]. This perturbation corresponds to the symmetry under field 
shifts of the evolution equation before the rescalings by k in eqs. (4) [39]. It is not 
expected to generate genuinely new physical behaviour, as it corresponds to a redundant 
operator [40] . The dominant non-trivial behaviour near the Ising fixed point is associated 
with the even perturbation. The two odd perturbations we discussed are related either 
to a shift of the source term (that is linear in a) or a field shift. All other perturbations 
have positive eigenvalues and are expected to become negligible in the infrared. 

In spite of the apparent absence of strong effects associated with the odd perturbations, 
the number of parameters that need to be fine-tuned in order to approach the fixed point 
is three in the absence of a Z 2 symmetry. In the model we are considering, with a bare 
potential given by eq. (2), we have three parameters at our disposal: m 2 , A, j. Moreover, 
the field value a* around which we expand the potential is not determined a priori If we 
fix it arbitrarily, we need to fine-tune all three of m 2 , A, j in order to approach the Ising 
fixed point. Equivalently, for an arbitrary value of j the fixed point can be approached by 
fine-tuning m 2 , A, a*. We have verified this conclusion through the numerical integration 
of the evolution equation. If the three-dimensional theory results from a more fundamental 
theory related to QCD, such as the quark-meson model, the parameters in the potential 
are a function of the temperature T and the chemical potential /i, while the source term 
j is proportional to the current quark mass. Then we conclude that for given j the fine- 
tuning of T and \x can lead to a second-order phase transition. By varying j a line of 
second-order phase transitions appears in the phase diagram. 

We are interested in the universal equation of state near the phase transition. For 
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this reason we need to discuss the evolution of the potential after the fixed point has 
been approached. The fine-tuning of the bare parameters guarantees that the system will 
remain near the fixed point for a long "time" t. However, if the fine-tuning is not perfect, 
the perturbations with negative eigenvalues will start growing eventually. We can solve the 
evolution equation (23) if we split Uk into even and odd parts as Uk = u% + u% and assume 
u t ^ u k- This assumption is a good approximation, especially in the vicinity of the fixed 
point, because the symmetric part includes the fixed-point contribution uj. The even part 
satisfies eq. (23), while for the odd part we use the ansatz u%(a) = f 2 (t) a + f 3 (t) u^(a). 
This leads to the equations 

f = -\*-<i>f* < 25 > 

f = -5< 1 + ">A < 26 > 

whose solution is f 2 = c 2 k~ h l 2 Z k X ^ 2 , f 3 = c 3 k~ x l 2 Z l J 2 , with Z k the wavefunction renor- 
malization and c 2 , c 3 arbitrary constants. This solution is consistent with our previous 
discussion of the perturbations near the fixed point (where r\ is constant), but also de- 
scribes the final evolution of the potential away from it. 
For k — > we obtain for the effective potential 

U(a) = U(a) - U° = U s {a) + c 2 a + c 3 U s '(a) ~U s (a + c 3 ) + c 2 a. (27) 

We know that U s (a) is related to the universal equation of state of the Ising model through 
eq. (12). This implies that we can parametrize the equation of state similarly to eq. (20) 

dfl 5m 2 
S i=to= F ^ 5m ") + + + et'fzM, * = | ff + e |i//» . ( 28 ) 

where 5m 2 , e — > and \a\ takes values in a range of a few |e|. The universal function 
fz 2 { x ) is specified by the Ising universality class [21, 37]. The critical exponents are 
P = 0.33, 5 = 4.8. For an effective theory resulting from QCD, the parameters 5m 2 , e are 
functions of the temperature and chemical potential, while 5j corresponds to deviations 
of the quark mass from a constant value. 

The appearance of the function F 2 (e, 5m 2 ) can be understood as follows: For a po- 
tential expanded around a fixed value a^, the parameters c 3 (related to e), 5m 2 (used 
in the parametrization of the symmetric part) are functions of the parameters m 2 , A of 
the bare theory (or T, ji for QCD). The parameter c 2 depends on j* as well, but in a 
trivial way: c 2 = — j* + G(m 2 ,\). The choice of j* does not affect the renormalization 
flow because the evolution equation involves second functional derivatives with respect to 
the field. As a result, a term linear in a does not change through the evolution and can 
be added directly to the effective potential. On the other hand, if the bare potential is 
not Z 2 -symmetric its first derivative at the origin changes during the evolution and can 
become non-zero, even if it was zero for the bare theory. Expressing m 2 , A in terms of 
5m 2 , e leads to the expression c 2 = — j* + F 2 (e,5m 2 ). The fine-tuning of j* amounts to 
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choosing it so that it cancels F 2 (0, 0). This leads to eq. (28), where F 2 (0,0) = and 5j 
measures deviations of the source from j*. The function F 2 (e, 5m 2 ) can be determined for 
a given bare theory similarly to the case of eq. (20). A simple analytical calculation is 
not feasible, as perturbation theory is not applicable away from the tricritical point. The 
evolution equation has to be integrated numerically while performing a triple fine-tuning 
of the initial conditions so that the Ising fixed point is approached. 

As expected the odd perturbations do not lead to the introduction of a new universal 
function. They simply generate a field shift. Our result is in agreement with ref. [41]. 
In that study the bare potential was assumed to have a very simple form (a polynomial 
of fourth degree). Our result provides the generalization to the case of an arbitrary bare 
potential. 

Our conclusions are not expected to be altered when higher terms in the derivative 
expansion are taken into account. The form of the odd perturbations will remain the 
same. The first perturbation corresponds to a term linear in the field. The introduction 
of such a term does not modify the exact flow equation for the effective action [19, 21] 
before the rescaling by k in eqs. (4), as this equation involves second functional derivatives 
with respect to the field. As a result the implications of such a term are independent of the 
truncation. Similarly, the second odd perturbation is related to a redundant operator that 
leads to field shifts. Again, its effects are expected to be independent of the truncation. 

In analogy with the crossover for the equation of state between the form characteristic 
of the tricritical point (eq. (18)) and that of the 0(4) universality class (eq. (12)) in 
the case with j = 0, we expect a different crossover as a function of j*. For j* — > the 
universal equation of state near the critical point is given by eq. (20), while for large it 
takes the form of eq. (28). A crossover is expected between the two forms as j* is varied 
continuously. 

7 Implications for QCD 

Because of its universal form, eq. (28) is expected to describe the critical behaviour of all 
physical systems whose phase diagram includes a critical point similar to that of our toy 
model. For two-flavour QCD the connection has been made more explicit in ref. [25]. It 
was shown that a potential whose leading terms are given by eq. (2) results, at energy 
scales below the temperature, from the integration of the fermionic contributions in a 
model of quarks coupled to mesons. The subsequent decoupling of the pions at lower 
scales leads to a critical theory with only one massless field, the a field. The universal 
equation of state is given by eq. (28), with 5m 2 , e proportional to deviations of the 
temperature 5T and the baryonic chemical potential 5/j, from certain values. The source 
term 5j is proportional to the deviation of the quark mass from a constant value. 

It is obvious that in experimental situations the critical point of QCD can be ap- 
proached only along the surface 5j = 0. The details of the experiment (center of mass 
energy, type of colliding nuclei) determine the effective temperature and chemical poten- 
tial. Information from various experiments can be used in order to approach the critical 
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point along a curve e = e(5m 2 ). An important question is whether the universal properties 
of the critical system (such as scaling parametrized by critical exponents) are observable. 
The solution of eq. (28) with Sj = 0, e = e(5m 2 ) determines the location of the vacuum. 
The function F2(e(Sm 2 ) , 8m 2 ) is a regular function around the point (0,0) and can be Tay- 
lor expanded. (An explicit example is given for i*i(e, 5m 2 ) below eq. (20).) As ^(0, 0) = 
the leading term is c5m 2 , with c some constant. (We have assumed e = e(5m 2 ).) The 
solution of eq. (28) is \a + e(5m 2 )\ = \c5m 2 /D\ 1 / 5 , where D = /(0). This solution emerges 
because the critical exponents 5 = 4.8, (3 = 0.33 satisfy (35 > 1. 

The "unrenormalized" mass term d 2 U /da 2 (equal to the inverse susceptibility) scales 
as \a + e(5m 2 )\ s ~ 1 ~ |5m 2 |^ -1 ^' 5 . This implies that the effective exponent 7 e // takes the 
value 7 e // = (5 — l)/5 = 0.79. This should be compared with the standard value 7 = 1.24 
in the Ising universality class [21, 37]. For the exponent v parametrizing the divergence 
of the correlation length we find u e ff = 0.40. This is a consequence of the scaling law 
v = 7/(2 — rj) that relates u, 7 and the small anomalous dimension rj = 0.036. The 
standard value of v in the Ising universality class is v — 0.63. 

It must be pointed out that the mass term may scale with the standard value for the 
exponent along certain paths that approach the critical point. For example, for e = 0, 
5j = F 2 (0 : 5m 2 ) the standard scaling is obtained. However, this approach to the critical 
point is unphysical as the quark mass cannot be altered. Another possible path has 5j = 
and e(5m 2 ) given by the solution of the equation F(e, 5m 2 ) = 0. Again, such a fine-tuned 
path is unlikely to be realized experimentally. These conclusions are in agreement with 
previous studies of the tricritical and critical points, in which the arguments were based 
on scaling relations [35]. 

The effective values for the exponents that we derived above are smaller than the 
standard ones in the Ising universality class by approximately 40%. This implies that 
the divergence of quantities such as the correlation length or the susceptibility along 
the experimentally accessible paths is much slower than the naive expectation. This 
conclusion is a consequence of the fact that the external source (the explicit symmetry 
breaking term) cannot be altered, as it is proportional to the quark mass. The smallness 
of the effective exponents implies that the universal behaviour is not easily accessible. The 
critical point must be approached very closely before the divergences of the correlation 
length or the susceptibility become apparent. 

As a final comment we mention that the inclusion of the strange quark is not expected 
to modify our conclusions. It affects primarily the location of the line of critical points, 
while the general structure of the phase diagram remains the same. The general conclusion 
is that there is a direction in the phase diagram that is not under experimental control. 
It is related to the explicit breaking effects of the chiral symmetry (SU(3)l x SU(3)r in 
the case of three flavours) that are associated with the masses of the quarks. As a result, 
the experimental approach to a critical point is expected to follow a path that displays 
the effective exponents we derived above. 
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